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IHTRODUCTIOS 


Aerodynamic  ground  teat  facilities  for  simulating 
re-entry  flows  involve  the  use  of  very  high  operating 
pressures  at  high  temperatures  (Ref.l).  Under  these  condi¬ 
tions  the  equation  of  state  departs  markedly  from,  the  perfect 
gas  lav  due  to  compressibility  effects  and  excitation  of 
internal  energy  modes.  The  associated  thermodynamic  proper¬ 
ties  are  hence  also  strongly  influenced.  Optimisation  stu¬ 
dies  of  ground  test  facilities  require  that  such  gas  imper¬ 
fections  are  known  with  a  reasonable  degree  of  accuracy. 

This  report  describes  an  experimental  program  to  make 
accurate  measurements  of  pressure,  volume  and  temperature  of 
nitrogen  and  air  at  pressures  up  to  3000  atmospheres  and 
temperatures  up  to  3500  K,  a  range  of  conditions  in  which 
the  compressibility  effects  are  most  dominant. 

The  behaviour  of  dense  gases  may,  in  principle, 
be  predicted  using  quantum-statistical  mechanics.  However, 
the  mathematical  difficulties  involved  when  molecular  inter¬ 
actions  must  be  taken  into  account  are  so  formidable  that 
solutions  have  only  been  found  using  highly  simplified  mole¬ 
cular  models  (Ref. 2).  Kany  attempts  have  been  made  to  obtain 
semi-empirical  equations  of  state  which  are  valid  for  any 
dense  gas  by  the  application  of  Van  der  Waals'  principle  of 
corresponding  states  (Refs.  3  and  4).  These  depend  on  having 
available  knowledge  of  the  equations  of  state  over  the  range 
required.  Once  the  equation  of  state  is  defined,  all  thermo¬ 
dynamic  properties  can  be  obtained  from  existing  relations 
(Ref. 4). 
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Up  to  present  experiments,  providing  the  knowledge 
on  which  the  equations  of  state  are  based,  have  been  confin¬ 
ed  to  the  studies  of  gases  at  high  pressures  and  ordinary- 
temperatures,  where  measurements  may-  be  made  under  steady 
state  conditions.  A  large  amount  of  data  is  available  in 
this  regime.  The  contrary  is  true  at  high  temperatures 
because  of  the  problems  encountered  in  maintaining  the 
structural  integrity  of  the  containing  vessel.  These  problems 
can  be  overcome  by  rapidly  heat'ing  and  compressing  a  gas 
sample  by  a  piston,  and  taking  the  desired  measurements 
under  transient  conditions. 

The  driver  section  of  the  VKI  Piston  Driven  Shock 
Tube  has  been  modified  to  generate  samples  of  dense  high 
temperature  gases  and  instrumentation  has  been  developed  to 
measure  pressure,  temperature  and  density  (Refs.  5,  6).  A 
sodium  spectral  line  reversal  technique  is  used  to  measure 
the  temperature  of  the  gas.  The  density  of  the  gas  is 
determined  indirectly  by  measuring  the  piston’s  position 
during  the  compression  process.  This  measurement  involves 
the  use  of  an  eddy  current  displacement  transducer.  A  quartz- 
crystal  piezo  electric  transducer  is  used  to  measure  the 
pressure  of  the  gas.  More  details  of  the  instrumentation 
are  given  within  and  in  Ref. 6. 

The  report  describes  experiments  carried  out  to 
generate  thermodynamic  Btate  data  of  nitrogen  and  air  for 
an  extended  range  of  conditions  than  reported  in  Ref. 6. 
Particular  attention  is  paid  to  new  observations  and  new 
techniques  developed  recently.  The  results  for  nitrogen  are 
then  compared  with  the  EnkenhuB  and  Culotta  (Ref. 4)  equation 
of  state  model  which  agrees  accurately  with  the  tables 
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generated  'by  Grabau  and  Brahinsky  (Ref. 7)  and  the  results 
for  air  are  compared  with  the  tables  generated  by  Grabau 
and  Brahinsky  (Ref. 8).  These  latter  tables  were  developed 
at  AEDC. 
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2.  EXPERIMENTAL  TECHKIQUE 

2 . 1  The  test  gas  compression  system 

The  system  for  compressing  the  gas  is  illustrated 
in  Pig. 1 .  A  piston  constructed  of  hardened  steel,  aluminium 
and  nylon,  weighing  12.6  kg  and  of  length  2J0  mm  is  free  to 
move  in  a  barrel  of  92.8  mm  diameter  and  1.98m  in  length. 

The  barrel  is  attached  to  a  reservoir,  of  diameter  270  mm 
and  length  2.3  m  at  one  end  and  is  closed  at  the  other  end. 
Initially  an  aluminium  diaphragm,  typically  of  thickness 
0.4  mm,  retains  the  piston  at  the  reservoir  end  of  the 
barrel.  The  gas  under  test  is  introduced  into  the  barrel] 
after  evacuating  and  purging  the  barrel  with  the  use  of  a 
vacuum  pump,  and  the  reservoir  is  charged  with  air  up  to  its 
operating  pressure  Pq  at  temperature  Tq .  The  piston  is 
released  and  driven  down  the  barrel  compressing  the  test  gas 
from  an  initial  pressure  and  Tj,^  to  a  maximum  pressure 

P4 1>  and  temperature  T4f.  Both  T0  and  Ti^  are  equal  to  the 
laboratory  temperature.  The  piston  is  instantaneously  at 
rest  when  the  peak  pressure  is  achieved.  After  the  first 
compression  stroke,  the  piston  reverses  its  motion  and 
returns  to  the  reservoir  end  of  the  barrel.  It  continues  to 
oscillate  in  the  barrel  until  brought  to  rest  by  friction. 
Values  of  Pifji  up  to  3000  kg/cm2  can  be  achieved  depending  on 
the  initial  pressure  ratio  across  the  piston.  This  limitation 
is  imposed  by  the  strength  of  the  barrel  and  its  end  fittings. 

The  pressure  (P4),  temperature  {T4)  and  density 
(P4)  variations  of  the  nitrogen  are  measured  during  the 
first  compression  stroke.  Typically,  the  initial  barrel 
pressure  measured  by  means  of  either  a  mercury  manometer  or 


8 


VO 


FlG  t-  General  assembly  of  compressor 


S9-9t-Ul-0Q3V 


AEOC-TR-76-65 


a  Wallace  and  Tiernan  gauge  is  between  0.1  to  10  kg/cm2; 
the  transit  time  of  the  piston  from  its  initial  rest  posi¬ 
tion  to  its  "peak  pressure"  rest  position  is  about  50  msec. 
The  total  measurement  time  of  interest  is  of  the  order  of 
1  msec.  During  this  time,  the  piston  is  between  10  and  1*0  mm 
from  the  closed  end  of  the  barrel. 


2 . 2  Pressure  instrumentation 

The  pressure  variation  of  the  test  gas  during 
the  stroke  is  measured  using  either  a  Type  6221  or  Type  6201 
Kistler  piezo-electric  pressure  transducer  and  a  Kistler 
568/M5  charge  amplifier  and  recorded  on  a  Tektronix  Type 
502A  Oscilloscope.  The  transducer  is  mounted  in  the  wall 
and  near  the  end  of  the  barrel,  as  depicted  in  Fig. 2. 

The  oscilloscope  trace  is  photographed  with  a  Polaroid 
camera.  The  calibration  and  assessment  of  accuracy  of  the 
pressure  measuring  system  is  described  in  Ref. 6. 


2.3  Temperature  instrumentation 

The  temperature  of  the  test  gas  is  determined  by 
using  the  sodium-line  reversal  technique  (Ref. 9)  illustrated 
in  Fig.3.  A  single  source  effective  double  beam  system, 
described  in  full  in  Refs.  5  and  6,  iB  used. 

In  the  early  tests,  a  Philips  D-8 ,  6  volt,  16-17 
amp  rated  tungsten  ribbon  lamp  is  used  as  a  background  light 
source  for  the  sodium  line  reversal  technique.  The  maximum 
temperature  achieved  by  this  source  is  approximately  2300  K. 
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In  order  to  extend  the  range  of  the  experiments  to  3500  X, 
a  Spindler  and  Hoyer  Carbon  Arc  Source  was  acquired  and 
used  instead  of  the  tungsten  ribbon  source.  The  Carbon 
Arc  Source  provides  a  non-varying  standard  temperature  of 
3300  X.  The  light  is  emitted  from  a  crater  in  the  positive 
electrode  caused  by  the  phase  change  from  solid  carbon  to 
gas  in  turn  caused  by  the  arc  across  the  electrodes.  The 
emissivity  of  the  centre  of  the  crater  is  assumed  to  be  that 
of  a  black  body  and  equal  to  1 .  A  detailed  discussion  of 
such  a  light  source  is  given  in  Refs.  10  and  11.  The  temper¬ 
ature  of  the  source  can  only  be  altered  by  applying  neutral 
filters  between  the  Bource  and  the  test  section; 

The  Spindler  and  Hoyer  lamp  iB  powered  by  the 
Institute's  220  V  DC  supply  which  has  a  very  small  residual 
ripple.  A  resistor  of  27*5  ohms  is  placed  in  series  with 
the  source  to  supply  it  with  its  rated  55  V  and  6  amps 
across  the  electrodes.  The  positive  carbon  electrode  ( 6  mm 
diameter  by  200  mm  length  type  NORIS  homogeneous)  is  posi¬ 
tioned  horizontally.  The  negative  carbon  electrode  (7  mm 
diameter  by  125  mm  length,  type  H0RIS  homogeneous)  is  posi¬ 
tioned  vertically  at  a  distance  from  the  positive  electrode 
varying  between  2-5  mm.  The  voltage  across  the  anode  and 
cathode  during  operation  varies  with  distance  between  them. 
In  order  to  allow  continuous  operation  then  an  electric 
motor  is  fitted  to  feed  the  electrodes  into  their  correct 
position  as  they  are  burnt.  The  rate  of  feed  is  controlled 
by  means  of  keeping  the  voltage  across  the  electrodes  to 
within  pre-det ermined  limits.  Because  of  the  method  used  the 
feed  of  electrodes  is  not  smooth.  Two  other  sources  of 
uncertainty  were  examined.  Firstly  the  27,5  ohm  resistor 
becomes  heated  due  to  the  high  current  passing  through  it , 


13 


AEDC-TR -75-65 


causing  resistance  and  hence  supply  voltage  changes.  This 
effect  is,  however ,  compensated  by  the  continuous  feed 
system.  Secondly  the  displacement  of  the  positive  electrode 
by  several  millimeters  during  feeding  causes  a  tendency  to 
place  the  system  in  an  unfocussed  condition. 

The  performance  of  the  carbon  arc  and  the  optical 
system  vas  examined  using  an  Evershed  and  Vignoles  optical 
pyrometer  in  a  similar  manner  to  that  described  in 
Bef.6.  When  focussed  directly  on  the  crater  of  the  positive 
electrode,  the  optical  pyrometer  gave  a  temperature  of 
3770  K  with  a  reading  scatter  of  ±40  K.  This  scatter  is 
attributed  to  subjectivity  errors  in  reading  the  pyrometer, 
and  to  the  uncertainties  caused  by  the  variation  of  the  gap 
size  between  the  electrodes.  This  calibrated  value  compares 
well  with  the  expected  temperature  of  a  carbon  arc  source 
detailed  in  Ref. 10  of  3800  K. 

The  temperature  measured  at  the  test  section 
position  after  the  light  had  passed  through  the  measuring 
system  optics  vas  3k6o  K,  with  a  similar  scatter  as  before. 
Although  some  optically  flat  glass  slabs  were  tested  for 
their  application  as  neutral  filters  to  vary  the  background 
temperature  to  the  SLR  system,  these,  however,  were  not  used. 
It  was  found  that  the  source  could  be  used  to  provide  an 
adequately  large  range  of  conditions  below  3500  K  without 
such  filters. 

The  overall  uncertainty  of  the  system  with  the 
carbon  arc  background  source  vas  assessed  to  be  il.jjE,  i.e. 
similar  to  that  incurred  using  the  tungsten  ribbon  back¬ 
ground  source  at  low  temperatures  {Ref. 6). 
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2.  1*  Density  measurement 

The  density  variation  of  the  compressed  nitrogen 
is  determined  from  accurate  measurements  of  the  piston's 
position  relative  to  the  barrel  during  the  final  stages  of 
the  compression  stroke. 

Two  complementary  systems  are  used.  A  small  lead 
pin,  mounted  in  the  end  wall  of  the  barrel,  is  crushed  by 
the  piston.  The  length  of  the  crushed  pin  gives  the  displace- 
ment  L  of  the  piston  from  the  end  wall  at  peak  pressure. 
Away  from  peak  pressure,  the  piston's  position  is  monitored 
at  intervals  of  one  millimeter  by  a  displacement  transducer 
mounted  in  the  side  of  the  barrel.  The  pin  and  transducer 
positions  are  shown  in  Fig. 2.  The  latter's  position  prevents 
the  transducer  from  being  subjected  to  pressure  levels 
above  30  kg/cm2.  The  transducer,  a  commercially  available 
Vibroraeter  TW6-100/A  eddy-current  contactless  displacement 
transducer,  monitors  accurately  machined  grooves  on  an 
aluminium  sleeve  attached  to  the  piston.  The  transducer  is 
excited  by  a  Vibrometer  type  100  TRI/A  carrier  amplifier. 

Knowing  the  piston  position  relative  to  the  trans¬ 
ducer  ,  the  displacement  (L)  and  volume  ( V i, )  of  the  compressed 
gas  may  be  determined  from  the  geometry  of  the  barrel  and 
piston.  The  number  of  moles  M,  of  gas  in  this  volume  is 
known  from  a  measurement  of  the  initial  volume  Vi,^,  pressure 
P4.,  and  temperature  Ti,^ .  The  density  of  the  gas  is  given  by 


P4 


Vi* 
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Pi* 


P“i  Y“i 
RT4.  V4 


(1) 


where  R  is  the  gas  constant  per  mole.  In  deriving  Eq.  1  ,  it 
is  assumed  that  no  gas  leaks  past  the  piston  during  the 
compression  stroke. 

An  assessment  of  the  uncertainty  of  the  measur¬ 
ement  of  density  is  given  in  Ref. 6. 


2*  5  Compression  cycle  measurement  control  system 


During  the  final  stages  of  the  piston  compression 
stroke,  the  pressure,  temperature  and  piston  position  signalB 
are  recorded  simultaneously.  The  output  from  the  position 
transducer  triggers  the  sweep  of  the  oscilloscopes.  The 
vertical  amplifier  of  the  Tektronix  Type  531A  oscilloscope 
used  to  record  the  position  measurement  operates  in  the 
DC  mode  and  is  adjusted  so  that  a  signal  is  recorded  only 
when  the  grooves  on  the  piston  are  traversing  the  transducer. 
The  time  scale  of  the  trace  is  extended  by  using  the  four 
alternate  trace  capability  of  the  oscilloscope's  type  M 
plug-in  unit.  A  four  trace  "raster"  display,  which  has  a 
calibrated  sweep-back  time  of  30ysec#  is  obtained.  The  oscil¬ 
loscope  monitoring  pressure  and  temperature  signals  operates 
in  the  single-sweep  mode  to  prevent  recording  compression 
strokes  other  than  the  first.  Base  lines  on  the  traces  are 
added  just  before  each  test.  To  correlate  the  time  scales  of 
the  two  oscilloscopes,  a  one-kHz  square  wave  of  small  ampli¬ 
tude  is  superimposed  upon  the  pressure  measurement  and  piston 
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position  measurement  signals . 


2 • 6  Summary  of  measurement  uncertainties 

The  following  table  summarises  the  results  of 
the  uncertainty  study  of  the  measurements  at  peak  conditions 
for  which  detailed  discussion  is  given  in  Ref. 6. 


Measurement 

Femarks 

Error 

Pressure 

Transducer  and  charge 
amplifier  errors  assumed 
small  since  directly 
calibrated . 

Oscilloscope  calibration 
and  trace  reading. 

M 1  .0!E 

Temperature 

Error  in  pyrometer  reading 
including  that  due  to  50  Hz 
ripple  and  feed  mechanism 
Trace  reading  1JS 

•\<±0 . 2% 

Position 

Lead  pin  measurement 

Timing 

MO,  3% 
*±1.0* 

Synchronis  at  ion 

- — — - - 

MO. 5% 

Overall  error  (as  a  sum  of 

above)  :  i* .  5 Jg 


In  addition  to  this  table,  uncertainties  arise  for  A>500 
(see  Section  U.3.1)  and  when  a  shock  is  present  (see  Section 
^.3.2). 
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Because,  smaller  deflections  away  from  the  peak 
conditions  are  measured,  and  furthermore  synchronisation  is 
more  difficult  in  these  regions,  the  sum  of  the  overall  errors 
is  expected  to  be  in  the  range  of  ±1%  to  10jE.  One  source  of 
uncertainty,  that  of  the  effects  of  impurities  on  the  experiments, 
remains  difficult  to  assess  as  discussed  in  Kef , 6 ,  More  will 
be  said  about  impurities  in  Section  4.5. 
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3.  COMPUTER  PROGRAM  AIDS  TO  THE  EXPERIMENTS. 


3 . 1  PiBton  compression  cycle  predictions 

Calculations  concerning  the  compression  of  the  gas 
have  been  carried  out  using  a  piston  cycle  program.  The  basic 
program  uses  a  4th  order  Runge-Kutta  numerical  method  to 
solve  the  system  of  two  ordinary  differential  equations 
resulting  from  Newton’s  law  describing  the  forces  on  the 
piston,  including  friction  during  the  cycle.  The  program  has 
been  modified  in  several  way.s  as  described  in  Ref.  6 
providing  aids  to  the  estimation  of  various  sources  of 
uncertainty  and  to  the  carrying  out  of  the  experiments. 

3 . 2  Data  reduction  programs 

Two  main  data  reduction  programs  have  been 
designed.  One  program  uses  inputs  of  the  measurements  of  the 
traces  as  read  directly  from  the  measuring  table  and  all  the 
carefully  determined  calibration  constants  to  calculate  the 
raw  data.  The  second  program  presents  the  results  to  allow 
comparison  with  appropriate  equations  of  state. 

1 .  Dat  a_reduct  £ 

The  values  of  the  state  parameters  are  generated  by 
introducing  into  the  program  the  co-ordinates  of  the  pressure, 
emission,  absorption  and  position  traces.  The  first  three  are 
introduced  at  equal  time  intervals j *the  latter  is  introduced 
as  the  time,  arbitrarily  referenced,  for  the  distance  trans¬ 
ducer  to  sense  the  equal  interval  distances  defined  by  the 
centre  of  the  teeth  and  the  grooves  of  the  piston. 
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The  folloving  information,  the  acquisition  of 
which  is  descnhed  in  Ref. 6,  is  also  provided  : 

-  oscilloscope  and  transducer  calibration  curves  and  values, 

-  lamp  voltage  and  current,  and  lamp  calibration  curve, 

-  transmission  coefficients  of  the  optical  system, 

-  crushed  lead  pin  length, 

-  corrections  for  barrel  expansion  and  piston  compression  due 
to  the  high  pressures  developed. 

The  reduced  data  iB  then  rearranged  in  equal 
interval  time  steps  by  using  a  four-point  interpolation 
procedure.  The  data  can  be  plotted  out  on  the  VKI  line 
printer  to  the  resolution  available  on  the  printer  itself. 

The  main  use  of  this  plot  was  to  give  a  quick  look  of  the 
data  (without  laborious  plotting)  and  to  indicate  where 
errors  in  transcribing  the  raw  data  from  the  traces  may 
have  occurred.  A  sophisticated  curve  plotter  for  presenting 
the  results  with  better  resolution  was  not  readily  available 
at  VKI. 


The  unitB  selected  to  present  the  reduced  gas  state 
data  were  kg/em2  for  pressure,  amagats  (defined  as  density 
of  the  gas,  i .  e .  M2  or  air,  at  standard  conditions*)  for  density 
and  K  for  temperature. 


onS  : 

Po  =  latm  ■  1.03323  kg/cm2  =  760  mm  Hg  =  1013214.6  N/m2 
Tq  =  273. 15  K 

Then  1  amagat  *  U U . 5 868  *  10~6  gm  -  mole/cm3 
i.e.  1  amagat  for  N2  ■  1. 21*903  x  10-3  gm/cm3 
and  1  amagat  for  air  =  1.291 1*3  *  10-3  gm/cm3 
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It  has  been  discussed  in  Ref.  6  that  despite  the 
electronic  synchronisation  of  traces  carried  out,  the  peak 
temperature  did  not  coincide  exactly  with  the  peak  pressure 
and  density.  Since  it  is  difficult  to  synchronise  the  data 
to  better  than  of  order  10  usee  it  was  suggested  that  the 
data  should  be  re-aligned  such  that  the  peaks  did  coincide. 
The  philosophy  for  this  is  that  provided  the  test  slug  is 
homogeneous,  it  is  difficult  to  conceive  that  the 
temperature  will  be  rising  whilst  the  pressure  and  density 
are  falling.  For  this  reason,  a  modification  was  made 

to  the  program  to  re-align  the  data.  This  was  carried  out 
by  finding  an  average  axis  of  symmetry  around  the  peak. 
ValueB  of  time  have  thus  been  re-calculated  with  respect  to 
this  common  axis  for  presentation  both  in  the  computer  plots 
and  in  the  second  part  of  the  data  reduction  outlined  below. 


2 .  Comparison_betveen_exgeri me ntal_data_and_ theoretical 


A  computer  program  was  devised  to  aid  the 

comparison  of  theoretically  determined  equations  of  state 

with  the  experimental  results.  It  calculates  for  each 

measured  point  (P,  p,  T)  the  following  parameters  : 

- *  T  )  i*e.  calculated  pressure  from  measured 
calc  seas  mess 

density  and  measured  temperature;  p  ,  (P  ,  T  ):  and 

calc  mess  mess 

T  -  (  p  „  „  „  •  .)■  Also  the  entropies  Sx/R  (p  ,  T  ), 

calc  mess  meas  1  mess  mess 

Sz/R(  „  »  ’I>_ _  )  »  S3/R(P  ,  p  )  are  calculated. 

E16&S  meas  meas  meas 


The  theoretical  thermodynamic  information  used  for 
nitrogen  was  the  Enkenhus-Culotta  equation  of  state  (Ref.  4) 
which  agrees  with  the  AEDC  tables  (Ref.  7)  to  within  a  few 
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per  cent  for  all  parameters  over  the  range  considered.  Two 
subroutines  give  pressure  and  dimensionless  entropy  as  a 
function  of  density  and  temperature.  For  air,  the  AEDC  tables 
(Ref.  6)  were  used.  A  four  point  interpolation  formula  gives 
pressure  and  dimensionless  entropy  as  a  function  of  density 
and  temperature.  In  each  case  a  I ewton-Raphson  searching  pro¬ 
cedure  is  used  to  calculate  P],  Tj ,  S2/  R  and  S3/R  when  the 
arguments  are  not  density  and  temperature. 

An  example  of  the  read  out  of  the  complete  data 
reduction  program  is  given  in  Table  1  placed  at  the  end  of 
the  report. 
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It.  RESULTS  AMD  DISCUSSION 


h . 1  Isentronicity  of  compression  cycle 

Further  observations  of  oscillations  around  the 
mean  pressure  variations,  similar  to  those  reported  of  Ref.  6 
have  been  sensed  for  lov  entropy  cases  as  illustrated  in 
Fig.  4.  These  cases  are  characterised  by  their  high  initial 
test  gas  pressure  requiring  greater  piston  velocities  to 
achieve  the  desired  conditions.  Further  calculations  using 
the  cycle  prediction  program  in  which  the  full  characteristics 
method  was  applied  to  the  flow  of  gas  in  front  of  the  piston 
(which  for  simplicity  was  assumed  perfect)  showed  that  weak 
shock  waves  can  be  developed  in  the  later  stages  of  the 
compression  for  these  high  piston  speed  cases  (i.e.  greater 
than  30  n/sec) . 

Further  examination  of  the  traces  indicates  that 
although  abrupt  pressure  steps,  indicative  of  shock  waves, 
are  not  seen  the  waves,  being  only  weak  ones,  are  indeed 
sensed  by  the  transducers  but  smoothed  by  the  limited 
response  (for  such  sharp  rise  times)  of  the  pressure  measuring 
system  .  Further  similarity  in  experiment  and  theory  was 
found,  in  that  these  weak  shocks  damp  out  quite  rapidly 
during  the  early  stages  of  the  expansion  process.  These 
shocks  could  be  reduced  in  strength  by  increasing  the  piston 
weight,  however  this  was  not  attempted  because  of  the 
limited  length  of  piston  that  could  be  introduced  into  the 
facility,  as  well  as  causing  a  problem  of  decreasing  the 
test  gas  sample.  Temperature  and  piston  position  measurements 
do  not  so  obviously  show  up  these  oscillations.  In  the  former 
casethis  is  because  the  oscillation  frequency  tends  to  be  of 


23 


AEDC-TR  -75-65 


the  same  order  as  the  chopping  system.  In  the  latter  case  the 
motion  of  the  piston  is  an  integral  function  of  the  pressure 
trace  and  hence  will  be  highly  smoothed. 

The  complete  calculations  have  not  as  yet  been 
carried  out  because  of  the  difficulty  in  introducing  the 
Fankine-Hugoniot  relations  to  account  for  flow  through  shock 
waves  which  are  curved  (since  they  are  moving  through  an 
accelerating  gas)  and  reflected  from  the  moving  boundary  of 
the  piston.  The  solution  is  possible,  but  considerable 
development  and  computing  time  is  required. 

The  main  consequences  of  this  formation  of  shocks 

are  : 

-  the  cycle  will  no  longer  be  isentropic;  small 
increases  in  entropy  will  occur; 

-  the  sample  will  be  non-homogeneous ,  causing 
larger  uncertainties  than  expected; 

-  transport  properties  ,  assessed  from  slopes  of 
data  from  expectedly  isentropic  compressions  and  expansions, 
will  become  inaccurate; 

-  analysis  of  the  basic  data  using  the  expectedly 
good  assumption  that  the  cycle  is  isentropic  should  be  treated 
with  caution. 

The  degree  to  which  the  formation  of  shocks  upsets 
the  sometimes  required  assumption  of  homogeneity  and 
isentropicity  has  not  as  yet  been  assessed.  However  further 
comments  are  given  in  Section  U.5* 
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^‘2  Some  observations  concerning  the  temperature 
measurement 


A  typical  temperature  measurement  trace  at  conditions 
of  high  temperature  and  pressure  is  illustrated  in  Pig. 5. 

It  is  seen  that  at  the  period  during  and  near  the  peak  temperature 
then  vi  is  equal  to  v3 .  As  pointed  out  in  Ref. 5  this  indicates 
that,  during  this  period,  maximum  absorption  of  the  light  from 
the  lamp  occurs.  The  equation  for  the  gas  temperature  : 


Tg  “  Tl  +  hv  ln  +  — 71 - >  (from  Eq. 1  in  Ref. 6) 

becomes 

!_  .  i_  .  k_  , 

Ta  tl  VI 

This  equation  may  be  shown  to  be  valid,  but  the 
measured  temperature  is  then  not  an  average  temperature  over 
the  optical  path  length  in  the  compression  chamber  but  an  average 
temperature  over  the  "optical  depth"  of  the  emitting  sodium  atoms. 
Provided  that  no  significant  gradients  exist  in  this  "depth", 
the  measured  temperature  is  still  indicative  of  the  test  gas 
temperature  and  no  additional  error  is  incurred  in  this  ease. 

This  condition,  v1  equal  to  v3 ,  occurs  at  combined  high  temper¬ 
ature  and  pressure  conditions  as  illustrated  in  Fig. 6.  The 
condition  has  an  advantage  of  making  the  temperature  measurement 
independent  of  the  lamp- iight  chopping  speed.  Advantage  has  been 
made  of  this  black-body’  method  to  measure  temperatures  in 
helium  for  temperatures  greater  than  2500  K  (Ref. 12). 
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P4j  a  0-399  kg/cm2 
T4-  a  294  K 
P0  a  15  kg/cm2 


FIG-5  -  Typical  high  temp* 


AEDC-TR -76-65 


NITROGEN 


to 

1_ 


1 

I _ 

Initial  pressure  (atm)  at  Te300  K 


0.1 

-j 


FIG- 6:  Regions  on  Motlfer  Chart  at  which  "Black  Body" 
measurement  occurs  -  a)  Nitrogen 
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Initial  pressure  (at  mi  at  300  K 


FIG-6:  Regions  on  Moilier  Chart  at  which  "Black  Body" 
measurement  occurs  *  bl  Air 
Black  Body  emission  is  thougth  to  extend  to 
the  dashed  line  if  no  chemical  reactions  or 

i 

Impurities  were  present. 
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If  maximum  absorption  does  not  occur  (i.e.  vi  f  v3 
near  the  peak)  then  the  traces  have  a  shape,  similar  to  that 
shown  in  Fig.  J,  In  all  the  nitrogen  traces,  good  symmetry 
of  the  emission  trace  is  obtained,  however  this  is  not  found 
in  air  cases.  In  all  cases  (nitrogen  and  air)  away  from  the 
peak  conditions,  asymmetry  of  the  absorption  trace,  however, 
is  obtained;  more  absorption  occurring  after  the  peak  than 
before.  This  asymmetry  has  been  discussed  in  detail  in 
Ref.  6,  in  which  high  level  turbulence,  impurities  arising 
from  ablation  etc.  during  the  expansion  part  of  the  cycle 
have  been  blamed  for  the  anomaly.  For  this  reason  more 
reliance  is  placed  by  the  authors  on  the  data  during  the 
compression  part  of  the  cycle. 


^ • 3 ■  Further  discussion  on  measurement 
uncertainties 

^•3.1  Uncertainties  arising  from  small  volumes 
of  gas  samples 

At  high  entropy  level  cases,  at  high  temperature 
conditions  the  minimum  volume  of  the  compressed  sample 
becomes  so  small  that  the  piston  approaches  the  end  of  the 
tube.  The  resolution  of  the  sensors  measuring  piston  motion 
then  approaches  the  value  to  be  measured,  hence  inaccuracies 
arise.  Inaccuracies  become  important  when  the  piston 
approaches  to  less  than  3  ram  from  the  end  wall  equivalent 
to  a  compression  ratio  *  of  approximately  500  for  N2  and  air. 


X 

The  compression  ratio  can  be  calculated  from  the  data 
reduction  program  as  the  density  measured  during  a  test 
to  the  initial  density. 


30 


AEDC-TR-7E-85 


P4j  *  0-399  kg/cm^ 
T4-  =  295  K 
P0  *  10  kg/cm^ 


NITROGEN 
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The  enclosure  for  the  temperature  measurements  has 
a  relatively  small  volume  i.e.  equivalent  to  a  length  of 
tube  of  approximately  0 . OU  cm.  The  limit  of  A  =  500  is  shown 
on  the  H-S  diagram  illustrated  in  Fig. 8.  This  limit  aligns 
approximately  with  the  T4  =  2700  K  line.  This  finding  tends 
to  nullify  the  usefulness  of  the  carbon  arc  light  source 
for  extending  the  range  of  temperature  to  3500  K  using  this 
particular  facility. 

^•3.2  Uncertainties  due  to  the  presence  of 
shock  waves 


As  discussed  in  Section  1* .  1  there  is  evidence  that 
shock  waves  are  formed  at  the  cases  in  which  the  piston  speed 
is  the  highest,  i.e.  low  entropy,  high  initial  pressure 
conditions.  In  the  data  reduction  the  philosophy  was  taken  to 
take  the  mean  reading  through  the  oscillations  occurring 
especially  in  the  pressure  traceB.  In  most  cases  the  oscilla¬ 
tions  are  no  larger  than  the  trace  width  and  hence  for  these 
cases  the  error  will  not  be  large.  The  largest  fluctuations 

-5?  &t  the  conditions  S/R  ■  22.5,  pressure  1660  kg/cm^, 
temperature  1660  K  and  219  amagats  (as  shown  in  Fig.b).  The 
regions  in  which  piston  speeds  are  greater  than  30m/see  correspond¬ 
ing  to  shock  formation  are  presented  in  Fig. 9.  In  nitrogen,  the 
limit  aligns  approximately  with  the  S/R  *  22.5  line  and  in  air  with 
the  S/R  =  23.5  line.  However,  the  overall  error  does  not  exceed 
5%  if  these  limits  are  decreased  to  the  S/R  ■  22  and  23  lines 
respectively  since  the  shock  is  very  weak  (less  than  0.5%)  for 
entropy  values  larger  than  22  and  23  respectively. 

Presentation  of  Results 


X 

Table  2  summarises  the  tests  in  which  measurements  of  the 
state  of  the  dense  high  temperature  gases  were  made.  The  initial 


x 

presented  at  the  end  of  the  report. 
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10  1  0.1 

L. _ 1 _ I 

Initial  pressure  (atm)  at  T&300  K 

PIG  B:  Conditions  on  MoLlier  Chart  at  which  X>  500 
a)  Nitrogen 
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'£ _ 1 _ 2« 

Initial  pressure(atm)  at  300  K 

FIG  8:  Conditions  on  Mollitr  Chart  at  which  A>  500 

b)  Air 
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S _ ! _ 5' 

Initial  pressure  (atm)  at  7=300  K 

FIG-9a- Region  where  a  weak  shock  forms  during 
the  compression  stroke  in  nitrogen 
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Initial  pressure(atm)  at  300  K 


FIG-Ob-Region  where  a  weak  shock  forms  during 
the  compression  stroke  in  air. 
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conditions  and  the  maximum  measured  values  of  pressure  , 
density  and  temperature  are  there  presented.  The  range  of  the 
tests  carried  out  in  this  series  i6  illustrated  in  an  H-S 
diagram  (Fig. 10)  by  plotting  the  range  of  pressures  achieved 
during  every  test  at  the  initial  dimensionless  entropy  of 
the  sample.  Because  avay  from  peak  conditions,  errors  are  greater 
than  4.5JS,  vertical  thick  bars  have  been  used  instead  of  lines. 

The  data  reduction  of  all  these  tests  using  the  method  described 
in  Section  3.2  and  plots  of  results  has  been  carried  out,  and 
copies  of  the  very  extensive  complete  set  of  results  rest  at 
VKI  and  AEDC. 

The  plots  of  peak  measured  values  of  pressure,  density 
and  temperature  against  their  values  calculated  from  the  other 
two  measured  parameters  {equivalent  to  Figs.  21,  22,  29,  30  of 
Bef.6)  revealed  only  that  even  larger  differences,  but  with 
similar  trends,  are  seen  over  this  extended  range  of  conditions. 

It  was  considered  that  no  usefulness  would  be  gained  by  plotting 
this  extensive  data  in  this  report,  in  view  of  more  useful  plots 
presented  in  the  next  section. 

4 . 5  Discussion 

The  aim  of  this  study  is  to  provide  new  experimental 
data  to  compare  with  proposed  models  of  the  equations  of  state 
at  conditions  well  outside  those  investigated  in  detail 
previously.  The  accuracy  of  these  theoretical  equations  of  state 
remains  unverified.  In  order  to  determine  the  state  of  a  gas, 
at  least  three  parameters  must  be  known.  In  .these  experiments 
pressure  and  temperature  were  measured  directly  and  density 
inferred  from  measurements  of  volume  of  a  known  constant  mass 
of  gas.  Great  attention  was  paid  to  keeping  the  uncertainties 
as  small  as  possible  ,  and  to  ensuring  the  quality  of  the 
sample.  When  all  of  the  information  has  been  obtained  severe  problems 
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Initial  pressure(atm)  at  300  K 


Ft  G - 10:  Summary  of  conditions  tested  -  b)  Air 
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exist  in  selecting  a  method  of  comparing  the  three  parameter 
experimental  results  with  the  models.  One  of  the  chief 
difficulties  is  to  ascertain  where  the  errors  lie  if  there 
are  discrepancies  arising;  they  may  he  in  the  development 
of  the  sample »  the  measurements  of  properties  of  the  sample 
or  in  the  actual  model  itself.  Another  difficulty  that  may 
arise  is  that  even  if  there  is  agreement  to  within  the 
assessed  uncertainty  of  the  experiment,  there  may  be  the 
chance  that  an  error  in  one  parameter  may  be  offset  by 
another.  However  in  the  latter  case,  such  a  state  of 

will  be  unlikely  if  agreement  is  found  in  many  tests 
taken  over  a  broad  range  of  conditions. 

The  approach  taken  in  the  early  stages  of  this 
research  (Ref.  6)  was  to  compare  the  measurement  of  one 
parameter  (either  P,  p  or  T) ,  with  the  calculation  of  it 
from  the  other  two  measured  parameters  using  the  model  for 
which  verification  is  desired.  The  problem  with  this  is  that 
if  differences  do  occur  then  it  is  difficult  to  decide 
whether  it  i b  due  to  an  error  in  the  model,  or  one  error  or 
an  accumulation  of  errors  in  the  experiment.  In  fact  as 
was  pointed  out  in  Ref.  6  systematic  deviations  from  the 
models  are  indeed  obtained  at  conditions  typically  above 
1700  K  and  1400  kg/em2  for  S/R  between  22  and  23.5.  The 
results  over  an  extended  range  (S/R  from  21.5  to  25.5)  show 
even  larger  deviations  as  noted  in  Section  4.4.  Air  tests 
show  larger  and  less  systematic  deviations  than  for  the 
nitrogen  tests. 

A  parameter  known  fairly  accurately,  but  as  yet 
unused  in  the  comparison  with  the  real  gas  equation  of  state 
is  the  entropy.  The  compression  is  essentially  isentropic,  such 
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that  the  initial  value  of  entropy,  which  can  he  defined  very 
accurately  from  the  initial  running  conditions,  will  remain 
constant  throughout  the  compression  and  expansion.  Small  rises 
in  entropy  are  expected  in  the  low  entropy  cases  during  a  test 
due  to  the  formation  of  weak  shocks  (see  Section  k.l)  and 
energy  losses  caused  by  heat  transfer.  In  order  to  use  this 
information  then  the  data  has  been  plotted  in  different  forms. 
Figs. 11  and  12  illustrate  for  all  entropy  level  cases  for 
nitrogen  and  air  respectively  measured  peak  temperature  plotted 
againBt  peak  pressure  and  compared  against  that  predicted  by 
the  representative  model  equation  of  state.  In  Figs.  13  and  Ik 
all  entropy  level  cases  are  plotted  for  peak  pressure  and  peak 
density.  In  Fig,  15,  peak  temperature  is  plotted  versus  peak 
density  for  two  entropy  cases.  It  is  seen  that  for  all  cases 
in  which  the  temperature  measurement  is  included  (illustrated 
by  Figs  - 1 1 »  12  and  15) »  then  deviations  away  from  the  equation  of 
Btate  models  are  produced.  It  should  be  noted  that  such  deviations 
were  seen  for  all  entropy  levels.  However,  it  is  seen  that  for 
all  cases  in  which  the  temperature  measurement  is  avoided  (given 
in  full  in  Figs.  13  and  Ik)  very  good  agreement  with  the  models 
are  produced,  except  for  high  entropy,  high  compression  ratio 
cases  when  large  uncertainties  are  expected  (see  Section  k.3.2). 

It  is  pointed  out  that  this  latter  plot  is  less  sensitive  to 
changes  in  entropy  than  in  the  previous  two  presentations  given. 
However  close  scrutiny  leads  to  the  conclusion  that  better 
correlation  is  indeed  obtained  in  this  plot. 

One  interpretation  of  these  observations  is  that 
the  temperature  measurement  is  in  error.  Possible  connected 
conclusions  arising  from  thiB  interpretation  and  from  the 
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FIG -l|  -  Plot  Of  peak  temperature  versus  peak  pressure 
In  n itrogen  • 

•Experiment,  —Equation  of  state  (Ref-7  J  at  given  S/R . 
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FIG- 12”  Plot  of  peak  temperature  versus  peak  pressure 
in  air- 

•  Experiment, —  Equation  of  state  (Ref- 8)  at  given  S/R. 
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FIG  -13- Plot  of  peak  density  versus  peak  pressure 
in  nitrogen. 

w  Experiment  ,  —  Equation  of  state  (Ref-7)  at  given  S/R- 
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1) 


FIG-14-  Plot  of  peak  density  versus  peak  pressure 
in  air- 

«  Experiment  , —  Equation  of  state  (Ref-8)  at  given  S/R- 
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FIG- 15*-  Plot  of  peak  temperature  versus  peak  density- 
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good  agreement  between  experiments  and  the  equation,  of  state 
models  are  that  : 

a)  the  density  and  pressure  measurements  are 

accurate . 

b)  The  compression  tube  produces  a  sample  of 
dense  high  temperature  gas  which  is  sufficiently  pure  and 
homogeneous  as  to  match  the  accuracy  of  the  instrumentation 
measuring  it. 

c)  The  equation  of  state  models  selected  are 

accurate . 


The  question  now  remains  as  to  where  lies  the 
error  in  temperature  measurement  in  this  hypothesis.  The 
lamp  and  optics  are  very  carefully  calibrated  and  in  theory 
there  appears  to  be  little  uncertainty  in  the  concept  of 
using  sodium  as  the  radiating  substance  in  these  high  density 
gas  samples.  Various  tare  tests  were  carried  out  to 
verify  that  no  errors  were  introduced  by  facility  vibrations. 
However,  no  direct  dynamic  calibration  of  the  method  under 
operating  conditions  is  available  at  present.  The  only  possible 
means  of  calibration  would  be  to  apply  another  temperature 
technique  at  the  same  time  as  the  sodium  line  reversal 
technique.  A  possible  reason  for  errors  in  temperature  is 
that  impurities  {not  occurring  in  large  enough  concentrations 
to  affect  the  overall  thermodynamic  characteristics  of  the 
gas  Bample)  may  be  affecting  the  sodium  emission  and  absorption 
of  the  light  from  the  lamp.  Little  is  known  about  the  effect 
of  impurities  on  the  SLR  method  under  conditions  of  high 
density.  To  the  authors'  knowledge  this  is  the  first  attempt 
made  and  the  effect  may  be  rather  sensitive  to  accuracy. 

For  all  nitrogen  tests,  and  the  lover  entropy 
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level  air  testa,  the  temperature  is  lover  than  expected 
(if  pressure  and  density  are  measured  accurately).  Equation  2 
illustrates  that  to  give  this  result,  either  the  light 
emission  is  too  lov  or  the  light  absorption  is  too  high. 

The  decrease  in  v3  (which  is  a  function  of  absorption)  to 
give  a  1  %  increase  in  temperature  is  found  to  be  50  %  at 
2200  K.  Similarly  the  decrease  in  vl  (which  is  a  function 
of  the  light  emission)  to  give  a  2  %  decrease  in  temperature 
is  found  to  be  50  %  at  2200  K. 

Several  characteristics  (noted  before)  may  provide 
some  doubt  about  the  method  of  measuring  temperature.  One 
is  that  the  absorption  traces  (and  sometimes  the  emission 
traces)  are  unBymmetrical  with  respect  to  peak  conditions 
in  which  more  absorption  occurs  during  the  expansion  than 
during  the  compression  at  equivalent  pressure  and  density 
conditions  (Section  k.2).  However,  in  view  of  the 
insensitivity  of  the  temperature  to  absorption  pointed  out 
in  the  last  paragraph,  this  is  thought  not  to  provoke  large 
errors.  The  other  concerns  the  observation  of  emission  in 
the  sodium  D-line  bandwidth  even  when  no  sodium  azide  was 
placed  in  the  tube  and  the  test  section  cleaned  thoroughly 
(See  Ref.  6).  Sophisticated  time  resolved  spectroscopic 
photographic  techniques,  beyond  the  capability  of  VKI ,  are 
required  to  understand  this  problem,  Even  after  such  a  study 
(likely  to  be  very  extensive)  it  may  not,  even  then,  be 
feasible  to  develop  the  sodium  line  reversal  technique  to 
sufficient  accuracy,  hence  it  is  thought  more  cost-effective 
to  examine  other  possibilities  for  measuring  temperature. 

Another  interpretation  of  the  poor  comparison 
of  experiments  with  the  equation  of  state  model  when  the 
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temperature  measurement  is  included,  but  good  comparison 
vhen  the  temperature  is  suppressed  may  be  again  associated 
vith  impurities.  In  this  case  the  impurities  may  be  so  large 
that  the  quality  of  the  gas  sample  is  impaired  but  the 
pressure,  temperature  and  density  measurements  are  accurate. 
Less  weight  is  given  to  this  interpretation  for  at  least 
the  nitrogen  tests  by  the  uniformity  of  the  correlations  of 
Fig.  13.  This  interpretation  has  more  weight  for  the  air 
tests  since  some  inconsistency  of  results  is  seen. 
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5  .  CONCLUSIONS 

The  compression  tube  has  been  shown  to  be  capable 
of  generating  samples  of  air  and  nitrogen  test  gases  up  to 
3000  kg/cm2  and  3500  K.  Results  of  measurements  of  pressure 
using  a  piezo-electric  sensor,  temperature  using  a  sodium  line 
reversal  technique  and  density  using  an  eddy  current  distance 
transducer,  show  discrepancies  on  comparison  with  advanced 
equation  of  state  models.  By  using  knowledge  of  the  entropy 
of  the  sample  it  was  illustrated  from  data  correlations 
that  the  discrepancies  most  probably  arise  from  errors  in 
temperature  measurement.  Such  errors  are  surmised  to  occur 
through  the  sensitivity  of  the  spectroscopic  technique  to 
slight  impurities  at  high  density  conditions.  Tentative 
conclusions  arising  from  these  successful  correlations  are 
that  the  transient  method  of  generating  dense  high  temperature 
samples  is  suitable  and  the  pressure  and  density  measurements 
and  the  gas  models  are  accurate.  Further  experiments  using 
a  non-spectroscopic  technique  for  measuring  temperature  are 
necessary,  however,  to  confirm  this  satisfactory  conclusion 
and  reject  other  possible  but  less  satisfactory  interpretations 
of  the  data. 


An  attempt  was  made  to  increase  the  temperature 
measuring  capability  to  3500  K,  but  in  view  of  the  preceding 
remarks,  and  the  fact  that  large  inaccuracies  in  the  density 
measurement  are  expected  at  such  high  temperatures  in  this 
facility,  the  usefulness  of  this  phase  of  the  study  was 
nullified. 
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TABLE  2  :  SUMMARY  OF  TESTS 


a)  nitrogen 


S/R-21.5  Nitrogen 


P 

max 

T 

nax 

P 

nax 

Test  No 

Kg/cm* 

K 

(Amagat } 

70 

2390 

1323 

73 

I960 

1390 

253 

75 

1770 

1 400 

235 

7  6 

2250 

1420 

265 

77 

2580 

‘  1420 

285 

S/R-22 

Nitrogen 

T8 

1220 

1450 

177 

79 

1215 

1320 

174 

60 

1900 

1420 

225 

81 

2560 

1435 

262 

82 

2720 

1625 

24? 

83 

3050 

1670 

256 

84 

2385 

1  44o 

245 

85 

2660 

1560 

242 

87 

3435 

1675 

301 

S/R-22 . 4 

Nitrogen 

14 

1955 

1795 

204 

15 

1168 

1567 

156 

20 

2394 

1830 

237 

26 

1200 

1555 

157 

S/R-22. 8 

Nitrogen 

5 

2240 

i960 

202 

6 

733 

1600 

109 

6 

1070 

1705 

136 

19 

2922 

2126 

234 

S/R=23. 1 

Nitrogen 

10 

1460 

1930 

154 

12 

2520 

2090 

206 

21 

3116 

2275 

187 

24 

2528 

2124 

203 

25 

2712 

2190 

212 
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8/R»23.5  Nitrogen 


Test  No 

P 

max 

Kg/cm2 

T 

max 

K 

^max 
(JLmagat ) 

13 

3020 

2320 

220 

23 

I960 

22 1*0 

163 

1 6U 

9b0 

1700 

100 

185 

1580 

1970 

128 

186 

2250 

2030 

166 

187 

2720 

2110 

185 

186 

3060 

2165 

186 

191 

b55 

ll*60 

69 

S/R=2b 

Nitrogen 

192 

305 

11*95 

b5 

193 

313 

1500 

1*8 

195 

131 

1305 

25 

196 

217 

1320 

3b 

197 

383 

1565 

5b 

198 

760 

1760 

86 

199 

1  U 1 0 

1790 

125 

200 

1320 

1955 

122 

201 

201  0 

2180 

163 

202 

2120 

2095 

157 

203 

2^35 

2175 

178 

S/R-2  U . 5 

Nitrogen 

20b 

220 

1570 

3b 

205 

283 

1550 

38 

206 

535 

1830 

60 

207 

1075 

2150 

95 

206 

1630 

22U0 

12b 

211 

1 1+2 

1380 

2b 

230 

1*790 

3060 

232 
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S/R=25.  Nitrogen 


Test  No 

P 

max 

Kg/cm* 

T 

max 

K 

p 

max 

( Amagat ) 

212 

680 

1945 

.60 

213 

745 

2080  . 

70 

214 

1045 

2280 

95 

215 

1535 

2385 

105 

216 

1236 

2270 

95 

217 

410 

1813 

44 

218 

217 

1625 

28 

220 

100 

1395 

16 

231 

4460 

3255 

175 

S/R=25.5 

N itrogen 

221 

235 

1745 

25 

222 

680 

1950 

54 

223 

1035 

2345 

71 

224 

965 

2330 

70 

233 

4150 

3450 

155 

Air 

S/R=22. 

Air 

228 

2235 

2030 

275 

S/R-22.5 

Air 

88 

740 

09 

993 

90 

1055 

91 

1185 

92 

966 

93 

1025 

94 

945 

95 

790 

96 

1665 

1570 

128 

1713 

163 

1745 

170 

1775 

181 

1760 

162 

1765 

_ 

1655 

157 

1680 

144 

1665 

218 
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Teat  No 

P 

max 

Kg/cm2 

T 

max 

K 

p 

max 

( Amagat ) 

97 

1 1*  1  5 

1920 

195 

98 

1280 

1750 

188 

99 

1765 

1900 

211 

100 

2365 

1990 

252 

102 

1990 

1905 

210 

103 

2990  ' 

1985 

287 

S/R=23. 

Air 

1  12 

467 

1500 

- 

113 

560 

1590 

108 

111* 

755 

1660 

130 

115 

1200 

1770 

165 

116 

1500 

1915 

185 

117 

2770 

2020 

245 

118 

11*20 

1935 

185 

122 

370 

1575 

74 

S/R=23 • 5 

Air 

123 

785 

1710 

120 

12U 

830 

1775 

- 

125 

855 

1750 

127 

127 

1620 

2010 

180 

128 

2685 

2045 

235 

129 

281*5 

2045 

- 

130 

3085 

2045 

245 

S/R=2U. 

Air 

38 

1500 

2050 

155 

39 

785 

1850 

105 

1*0 

727 

1830 

100 

1*1 

955 

2130 

120 

50 

1290 

2285 

1  40 

51 

2510 

2120 

191 

56 

1690 

2180 

165 

132 

2290 

2040 

187 

133 

3230 

2020 

232 

135 

3  55 

1605 

65 

136 

340 

1620 

64 
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Test  Ho 


P 

max 

Kg/cm2 


S/R=24. 5 

Air 

42 

420 

44 

680 

45 

1060 

46 

1430 

46 

985 

U7 

100 

153 

210 

154 

325 

155 

510 

156 

3335 

S/R=25 • 

Air 

157 

155 

158 

200 

159 

425 

160 

535 

161 

540 

S/R=25 . 5  Air 


T 

max 

K 


1550 

1710 

1630 

T680 

■1420 

1620 

1410 


1870 
1980 
2115 
2250 
2260 
1  314-5 
1  44o 
1810 
1915 
2235 


1580 

168O 

1760 

1940 

2130 


■  p 
■  max 

( Amagat } 


89 

115 

74 


59 

45 


63 
86 
113 
12  4 
102 
23 
39 
50 
71 
208 


27 

32 

53 

62 

64 


163 

330 

164 

44o 

I65 

645 

166 

735 

167 

1065 

168 

925 

169 

920 

177 

285 

178 

120 

1720  4l 
i860  50 
2040  71 
2045'  68 
2200  95 
2130  67 
2180  67 
1780  35 
1600  18 
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P 

T 

T) 

Test  No 

max 

max 

max 

Kg/ cm2 

K 

( Amagat ) 

179 

123 

1720 

19 

235 

1+205 

3275 

250 

S/R«=26. 

Air 

170 

1+1*5 

2055 

36 

173 

61*5 

2090 

1*3 

176 

1*33 

1950 

36 

236 

1*875 

3265 

165 

s/r=26.5 

Air 

18 1 

6B5 

2115 

36 

182 

305 

2010 

21* 

237 

21*15 

3290 

81 
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LIST  OF  SYMBOLS 

Planck  constant 
enthalpy 

Boltzmann  constant 

displacement  of  the  piston 

number  of  moles 

pressure 

gas  constant 

entropy 

sodium  line  reversal 

temperature 

volt  age 

volume 

frequency 
density 
di ameter 

Subscripts 

reservoir 

calculated  value  of  entropy  using  p  and  T 
emission 

i 

calculated  value  of  entropy  using  P  and  T 
lamp  emission 

calculated  value  of  entropy  using  P  and  p 
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final  conditions 


initial  conditions 
lamp 

sodium  atoms 


